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Abstract 

We show that quantum effects can stabilize a soliton in a model with no 
soliton at the classical level. The model has a scalar field chirally coupled to 
a fermion in 1+1 dimensions. We use a formalism that allows us to calculate 
the exact one loop fermion contribution to the effective energy for a spatially 
varying scalar background. This energy includes the contribution from coun- 
terterms fixed in the perturbative sector of the theory. The resulting energy is 
therefore finite and unambiguous. A variational search then yields a fermion 
number one configuration whose energy is below that of a single free fermion. 

I. INTRODUCTION 

Scalar field theories can contain spatially varying (but time independent) configurations 
that are local minima of the classical energy. These solitons are found as solutions to the 
non-linear classical equations of motion. Sometimes a topological conservation law can be 
used to show that the soliton is absolutely stable because it is the lowest energy configuration 
with a given value of a conserved topological charge. When quantum effects are taken into 
account, the classical description must be re-examined. Now the spatially varying soliton 
configuration should minimize the "effective energy" which takes into account classical and 
quantum effect sQ. Since the effective energy for general configurations is difficult to compute, 
quantum effects are typically computed as approximate corrections to the classical soliton. 
In a non-renormalizable theory, these corrections are cutoff dependent and the model must 
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be redefined to include the cutoff prescription. The hope is that the energy of the soliton is 
slightly altered by quantum effects but its qualitative features remain. 

In this Letter we give an example of a quantum soliton that is not present in the classical 
theory alone. We examine a renormalizable model in 1 + 1 dimensions where a scalar field 
is Yukawa coupled to a fermion. Fermion number is conserved. The classical energy is 
minimized when the scalar field has a constant value v. There are no classical solitons. 
The fermion gets a mass m = Gv through the Yukawa coupling. We calculate exactly the 
fermion's properly renormalized one loop contribution to the scalar field effective energy. 
By "exactly" we mean to all orders in the derivative expansion, which is crucial since we 
consider configurations varying on the scale 1/m. We then show that for certain choices 
of model parameters — in particular with G large — we can exhibit a field configuration 
that carries fermion number and has energy below m. It cannot decay by emitting a free 
fermion. We search for the lowest energy configuration carrying fermion number using a few 
parameter variational ansatz. The soliton, which is the actual lowest energy configuration 
with fermion number one, is presumably not far from our variational minimum and has 
strictly lower energy. The soliton therefore has energy less than m and is absolutely stable. 

The idea that a heavy fermion can create as soliton is not new and has been explored 
previously |I10] . What we are adding to the discussion is the ability to exactly calculate the 
renormalized fermionic one loop effective energy for any spatially varying meson background, 
which is essential for demonstrating stability at the quantum level. Since we are working in a 
renormalizable theory, the counterterms in the Lagrangian cancel the cutoff dependent part 
of the sum over zero-point energies in the explicit evaluation of the effective energy, leaving 
a finite result. This result is unambiguous because we are able to fix the counterterms in 
the perturbative sector of the theory. Furthermore, we can choose model parameters to 
justify neglecting the one loop boson contributions and all higher loop contributions. Thus 
we conclude that in this 1 + 1 dimensional model, a heavy fermion can create a soliton. 



II. THE MODEL 

The model we consider has a two-component meson field <fi = (0i, <p2) coupled equally 
to Np fermions. We suppress the fermion flavor label but will keep track of the factor Np 
as necessary. The Lagrangian is C = Lb + Cp with 

C B = \d^-d^-V($). (1) 

where 

= I 

and 

C F = ~(i[*,Ftf] -G 0i + i[*,75*] <h)) • (3) 

(The reason for the commutators in eq. (|3[) will be explained later.) Note that with a set 
to zero, the theory has a global U(l) invariance 
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0i + ifo -> e iv (0i + ifc) and * -> e -^ 7s/2 ^ 



(4) 



Naively, one would imagine that spontaneous symmetry breaking occurs with a = 0. Then 
we could pick a classical vacuum, say c i = (v, 0), and expand the theory about this point. 
But in 1 + 1 dimensions, the massless mode that corresponds to motion along the chiral 
circle, ■ = v 2 , gives rise to infra-red singularities and there is no spontaneous symmetry 
breaking [[J. By introducing a ^ we have tilted the potential to eliminate the massless 
mode. For a large enough it is legitimate to expand about <fi c \. There are two massive 
bosons, which we call o and n, with m 2 = (A + a) v 2 and m 2 = av 2 . The classical bosonic 
theory governed by Lb has no classical soliton. 

The fermions get mass through their Yukawa coupling to 0. In the perturbative vacuum, 
expanding about c i, the fermion has mass m = Gv. One could imagine that various 
distortions of would affect the fermion spectrum. For example, one could keep 02 = and 
let 0i — > 0i(x) with lim^^ioo 0i (x) = v, but <pi(x) < v over some region in x of order w. 
Alternatively, one could keep • = v 2 , but let = v(cos8(x), sin6'(x)), where 9(x) — > as 
x — > — oo and 6(x) — ► 2ir as x —>■ +oo. Again the deviation of from c i occurs in a region 
of order w. In both cases, if w is of order 1/m, then there are bound state solutions of the 
single-particle Dirac equation associated with eq. (|3]) that have binding energies of order m, 
so that a fermion bound to the field has an energy below m. (Because of its topological 
properties, the latter configuration is especially efficient at binding a fermion [|J.) However, 
there is an energy cost from the gradient and potential terms. Still, considering just the 
single bound fermion and the classical scalar energy, we might expect a total energy below 
m for G large enough. 

Of course, \l/ describes a quantum field and any distortion of the background <fi{x) away 
from c i will cause shifts in the zero-point energies of the fermion fluctuations. To form a 
self-consistent approximation, we must compute the effect of these shifts as well, since they 
are of the same order in h as the bound state contribution. In general the sum over zero-point 
energies diverges. In order to proceed we must regularize and renormalize the calculation. 
We are working in a renormalizable field theory so we know that the counterterms that 
are implicit in Cb will cancel these divergences. We want to compare the energy of non- 
trivial configurations with the perturbative spectrum of the model, therefore we fix the 
counterterms by standard renormalization conditions on the Green's functions. The Green's 
functions are evaluated perturbatively so the counterterms have an expansion in Feynman 
diagrams. 

Regularization and renormalization of the sum over zero-point energies has been prob- 
lematic in the past. We work in the continuum, where the sum is replaced by an integral 
over scattering phase shifts. In Ref. we show that it is possible to analytically continue 
this integral to d spatial dimensions, where it converges. Then we are able to identify po- 
tential divergences with low orders in the Born expansion for the phase shifts, and, in turn, 
with specific Feynman diagrams. We subtract the low order Born terms from the integral, 
which then remains finite when analytically continued back to d — 1. We then add back in 
the corresponding Feynman diagrams, which combine with the counterterms in the usual 
way to yield a finite and unambiguous result in d = 1. In the next section we show how we 
evaluate this contribution to the energy of a static configuration. 
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III. THE ONE FERMION LOOP EFFECTIVE ENERGY 




We have written eq. (|3]) as a commutator to ensure that the Lagrangian is invariant 
under the charge conjugation operation ^ — > C^* and ((/>i,</> 2 ) — > (<f>i,— 4>2)- As a result, 
the vacuum energy gets contributions from both the positive and negative energy eigenvalues 
of the single particle Dirac Hamiltonian 

H[ ( f]=ia 1 ^ + G ((7 2 0! + az<p 2 ) . (5) 

Here we are using a Majorana representation of the Dirac matrices, 7 = 02, 7 1 = icrs, and 
75 = (Ti, which implies that C = 1. For one flavor of fermions the energy of the lowest 
energy state is 

£vac = - \ - u n + y^u n \ = -^^2 i^ n ' ' ( 6 ) 

n 

where the o/s are the eigenvalues of H[4>]. (For a charge conjugation invariant background, 
for each eigenvalue u there is an eigenvalue —u, so the two sums in eq. @ are the same, and 
E VSLC reduces to the sum over the Dirac sea, _E vac = Yjw <o oj n .) We will restrict our attention 

to background fields that obey <j)\(x) = <f>i(—x) and fafe) = — 2 (— x). In this case, H[<j)\ 
commutes with the parity operator P = <7 2 n, where II is the coordinate reflection operator 
that transforms x to —x. We can thus decompose the solutions of eq. (|E|) into separate 
parity channels. 

For a given background, we wish to evaluate eq. (H) and subtract from it what we get in 
the free case, <fi = (v,0). Following Ref. 0, we use the relationship between the change in 
the density of states and the derivative of the phase shift 

to write the change in the vacuum energy as 



where 



5 F (k) = 5+(u(k)) + 6 + {-uj{k)) + S-(u{k)) + 6-{-u(k)) . (9) 



Here u(k) = \Jk 2 + m 2 , 8± is the scattering phase shift for the positive (negative) parity 
channel, the {u>i} are the discrete bound state energy levels, and E ct is the counterterm 
contribution, which is fixed by renormalization conditions discussed below. The extra m/2 
reflects an important subtlety in one dimension: in the non-interacting case (5p(k) = 0) 
there are bound states exactly at the continuum thresholds, uj = ±m, which count as 1/2 
in the sum in eq. @ 0. Levinson's theorem, 

S±{m)+6±{-m) = Tr{n ± -±), (10) 
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relates the phase shift at threshold to the number of bound states n±, with threshold bound 
states again counting as 1/2. It allows us to rewrite eq. @ as 

1 f°° dk d 

AE F = "2 E (M " m ) " I 2^ ~ m) dk 6F{k) + Ect ' (U) 

which is convenient because it makes it clear that as we increase the background and a 
bound state appears, there are no discontinuities in AE F . 

Of course, AE F given by eq. (|TTD is formally infinite. For large k, the phase shifts go 
to zero like 1/k so the integral is divergent. To regulate this divergence and allow us to 
identify it unambiguously with specific Feynman diagrams, we have extended the method 
of dimensional regularization to the density of states written in terms of phase shifts. The 
details are presented in Ref. ||. Once continued to <i-dimensions, where all quantities are 
finite, we can identify the leading large k behavior of <5f(/c) with the contribution of the 
first Born approximation plus the piece of the second Born approximation related to it by 
chiral symmetry, which we call 5(k). We also identify it unambiguously with the coefficient 
of the Lagrangian counterterm, v 2 — <p ■ <f), evaluated by standard Feynman perturbation 
theory. The renormalization conditions that fix the counterterm in perturbation theory here 
translate into the statement that in evaluating eq. ( |TT| ) we should subtract 5(k) from Sp{k). 
After this subtraction the integral can be analytically continued back to d = 1 to give a 
result that is finite and unambiguous, 

I /* OO 7 7 7 

AE F = --£ (N -m)-J - (u(k) - m) - (s F (k) - 6(kj) , (12) 



where 



S(k) = ^r J dx(v 2 - 0(x) 2 J . (13) 

We can solve numerically for the phase shift Sp(k) for any background <p(x). We make 
use of the fact that P commutes with H[<f)] so that positive and negative parity channels 
are decoupled. The positive (negative) parity states obey 

ip±(-x) = ±a 2 ^±(x) (14) 

and therefore 

^ ±(0)OC (±i) • (15) 

Any state defined for x > and obeying one of the boundary conditions in eq. (|1"5| ) can be 
extended via eq. (|T4]) to the whole line, so we need only consider the half line x > with the 
boundary conditions eq. ([15]). Consider the free case, = (v, 0). For each energy u, both 
positive and negative, the right moving eigenstate of eq. (||) is 



UJ 



5 



and the left moving free eigenstate is 



where k = \Juj 2 — m 2 is positive. For backgrounds 4>{x) that are not everywhere equal to 

— * — * — * 

c i, we still impose the requirement that <f>(x) approaches (j) c \ as x gets large. For these 
non-trivial backgrounds we call fk(x) the eigenstate of eq. (f|) that approaches fl(x) as 
x — > oo and (p^{x) the eigenstate that approaches <p^_ k (x) as x — > oo. Note that tp° k1 </j° fc , 
and <y9_fc are not eigenstates of P. For x > let 

V±(x) = ± e 2 ^Vfc(z) (18) 

be the parity eigenstates of H[(j>] with energy to. This defines the phase shifts 5±(u). The 
27r ambiguity in this definition is resolved by requiring that the phase shifts be smooth and 
go to zero as lj — > ±00. Note that in the free case 5±(lj) = 0. For any value of lo, we can 
solve numerically for the eigenstates of eq. @ in both parity channels. Using eq. (jlSf) we 
can then extract the phase shifts. Our method for computing the phase shift actually allows 
us to resolve the 2n ambiguity for each u individually and has other numerical advantages, 
which are elaborated in Ref. H. 



For any value of k we can obtain 5p(k), so we can compute the integral in eq. fljj). To find 
the bound state energies we solve the eigenvalue problem numerically. Levinson's theorem 
tells us how many bound states to search for. For a fixed background 0, the numerical 
evaluation of eq. (|T2| ) can be done quickly and with high accuracy. This allows us to search 
over a class of 0's for the configuration with the lowest total energy. 



IV. THE TOTAL ENERGY 

We are interested in calculating the total one loop effective energy of a static configuration 
<f>(x). We take 4>(x) to be specified by a short list of parameters {0}. We measure energy 
in units of the fermion mass m = Gv and use a dimensionless distance £ = mx. In 1 + 1 
dimensions <p(x) and v are dimensionless. We rescale <p(x) by v so that <j){x) —>■ (1,0) as 
|£| — > 00, and define dimensionless couplings 

a ~ A 

a = — and A = — . (19) 

By this rescaling, using eq. ([I]) and eq. (0), we have 




m 

= v 2 £ cl (a,X,{Ci}), (20) 

where prime denotes differentiation with respect to £. 

The fermion one loop contribution to the energy arises from eq. (|J), which with mea- 
sured in units of v is 
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H[$] = m Ua^ + + a 3 cf> 2 (0) • (21) 

We see that a single fermion makes a contribution proportional to m and dependent on the 
variational parameters {Q}- This means that eq. (0) can be expressed as m£ F ({Q}). For 
Np flavors the one loop contribution is therefore Npm£ F {{0})- 

The boson one loop contribution comes from summing the square roots of the eigenvalues 
of the operator — + J^.q^. ■ Rescaling as before we find that the boson one loop energy can 
be written as m£ B ({Ci}) • Putting together the classical energy and the one loop energies 
we get 

= ^g cl (a,~X,{d}) +£ F ({Q}) + -rj - £ B (<5? A, {Ci}) + higher loops. (22) 
Npm Np Np 

For Np large we can neglect the boson one loop contribution relative to the fermion one 
loop contribution. Furthermore it can be shown that 1/v 2 counts boson loops. Taking both 
Np and v 2 large with the ratio fixed, we can neglect the higher loops entirely and all but 
the single fermion loop in eq. (0). Therefore we need only consider the contributions from 
£ c \ and £ F . 



V. THE FERMION NUMBER 

A non-trivial background <p(x) distorts the energy levels of the Dirac Hamiltonian eq. ©, 
possibly introducing single particle bound states (with positive and negative energy). We 
identify the lowest energy state of the system, the one with the all the negative energy levels 
filled, as the vacuum. If a level crosses zero as we locally interpolate between <p c \(x) and 
4>(x), this state will have non-zero fermion number. In particular, if <fi(x) circles <fi = (0,0) 
as goes from (1,0) at x = — oo to (1,0) at x — oo, then the vacuum state will carry 
non-zero fermion number provided that the scale over which varies, w, is much larger 
than the fermion Compton wavelength 1/m |7|,|2[]. In Ref. || we derive a formula for the 
fermion number of the vacuum, <2 vac , in terms of the positive energy phase shifts at k = 
and the number of positive energy bound states, 

Qvac = N f(1 l 5 + ( m ) + 5 ~ ( m )] +\~ n ->o) (23) 

where n^ >0 is the number of bound states with positive energy. The configurations we look 
at loop at most once around = 0, so <2 vac is either or Np. We are interested in states 
with fermion number Np. If Q vac = Np, then the state we want is the vacuum. If Q vac = 0, 
then we build the lowest energy state with fermion number iVV by filling the lowest positive 
energy level of eq. (|5|) with iVV fermions. Therefore, if Q vac = 0, £ F appearing in eq. (|22|) 
must be augmented by uj\ where mio\ is the smallest positive eigenvalue of eq. (|2ll) . 
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VI. RESULTS 



We want to look for background configurations <\> that can produce states with fermion 
number Np, and whose total energy is below Npm. From eq. ( |2"2] ) with S B neglected, we 
define 

B = £-£d + £ F -l, (24) 

IV p 

which is the energy of the fermionic configuration minus the energy of Np free fermions in 
units of mNp. For our numerical computations, we take the ansatz 

(j) x + i<p 2 = 1 - R + R e ie where 6 = tt(1 + tanh(f /w)) . (25) 

The two variational parameters are R and w. As £ goes from — oo to oo, <fi moves in a circle 
of radius R in the (<fi\, fa) plane, starting and ending at (1,0). The scale over which (f> 
varies is w. 

For fixed at and A we vary R and w until we produce the configuration with the smallest 
B. The results are shown in Fig. |l[ We see that it is possible to find a configuration whose 
total energy is below Npm. Since we are minimizing B subject to the constraint that is 
of the form eq. (p5]), we know that the true minimum of B in the fermion number Np sector 
also has an energy below Npm. This is the stable soliton. 

In Fig. we show the width, to so i, and the radius, _R so i, for the minimum energy configu- 
ration as a function of a for various values of A. Note that the size of the soliton grows like 
1/ 1 \f& as a goes to zero. In that region, i? so i approaches 1, so the <fi configuration approaches 
the chiral circle. In fact, the energy of the fermion number iVV soliton goes to zero as a goes 
to zero. However for a very small the bosonic quantum fluctuations restore the classically 
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FIG. 2. The width u; so i (left panel) and the radius i? so i (right panel) of the configurations (25) 
that minimize the total energy as a function of the explicit symmetry breaking a. Several values of the 
Higgs coupling constant A are considered and v/^/Np = 0.375. 



broken symmetry. Thus we can not trust our results for a very small and we do not believe 
that this large and light soliton is a reliable consequence of this model. For moderate values 
of a, where the width of the soliton is not controlled by 1/y/di, we do trust our results. 
For the value of v/y/Np shown in Fig. |2|, the model becomes trustworthy for a ~ 0.3. For 
further discussion of this point, see Ref. ||. 

We have developed and applied a variational technique for renormalizable quantum field 
theories through one loop order. Because we have applied unambiguous, standard perturba- 
tive renormalization procedures, we have been able to hold the theory (i.e. the renormalized 
masses and coupling constants) fixed, while searching over a variational ansatz. Here we 
have used these methods to demonstrate the existence of a stable fermionic soliton, stabilized 
by quantum effects, in a model with no soliton at the classical level. The result suggests 
that similar phenomena might persist in 3+1 dimensions, and no obstacles stand in the way 
of generalizing the method to that case. 
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